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INTRODUCTION 


Projectiles  are  traditionally  analyzed  and  designed  by  using  theory  of 
elasticity  ,  H-r2-r3-} .  In  this  approach,  the  entire  projectile  body  is 
considered  to  be  loaded  within  the  elastic  region.  However,  under  actual 
firing  conditions,  the  equivalent  stress  in  some  region  of  the  projectile 
is  much  greater  than  the  yield  strength  of  the  projectile  material.  Hence, 
plastic  flow  will  be  encountered  in  the  projectile  body.  The  present  trend 
in  stress  analysis  of  weapon  and  ammunition  components  is  to  take  into 
account,  in  a  more  rigorous  fashion,  the  complex  phenomena  of  plastic 
flow.  This  is  due  to  the  necessity  of  designing  for  maximum  stress  to 
achieve  an  optimal  design.  In  this  investigation,  a  more  rigorous 
nonlinear  technique  is  developed  in  order  to  predict  the  Inelastic 
deformation  and  stress  distributions  of  a  projectile  subjected  to  actual 
firing  conditions.  Both  nonlinear  material  response  and  geometric  nonlin¬ 
earity  have  been  taken  into  consideration.  Nonlinearity  of  material  prop¬ 
erties  has  been  taken  into  account  by  use  of  theories  of  plasticity. 
Geometric  nonlinearity  has  been  considered  by  use  of  the  finite  element 
approach.  A  complete  inelastic  stress  analysis  of  a  30  mm  XM  TP  projectile 
has  been  conducted.  The  intent  of  this  investigation  is  to  identify 
potential  design  flaws  and  critical  regions  of  the  projectile  under  actual 
firing  environment. load  history  was  defined  which  subjected  the  pro¬ 
jectile  to  the  loads'at  various  points  on  the  barrel  travel  curve.  A  peak 
pressure  of  60,000  psi,  peak  acceleration  of  4.8  x  I07  in. /sec2  and  a  maxi¬ 
mum  spin  of  6,355  rad/sec  were  considered. 
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CONSTITUTIVE  EQUATIONS  j 

;  I 

The  primary  objective  of  this  Investigation  Is  to  develop  a  more 
rational,  nonlinear  elastoplastic  method  for  analyzing  and  designing  a 
projectile  to  achieve  in-bore  structural  integrity.  The  incremental 
stress-strain  relations  associated  with  the  Von  Mises  yield  criterion 
N  obeying  the  kinematic-hardening  law  for  work-hardening  materials  will  be 
used.  The  kinematic  hardening  law  in  the  incremental  theory  of  plasticity 
was  originally  proposed  by  Prager  (4],  for  the  deformation  of  isothermal 
►  solids •  Later,  Prager  [5]  extended  these  formulations  to  the  nonlsothermal 
i>  condition  for  rigid  work-hardening  solids.  Recently,  Chu  [6,7],  expanded 
Prager' s  work  for  solids  of  more  general  deformation  state.  A  brief  out¬ 
line  of  the  constitutive  equations  used  in  the  analysis  is  now  given.  ' 
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The  total  increment  strain  tensor,  dc^«  is  assumed  to  be  the  sum  of 
elastic  part,  de®,  ,  plastic  part,  dejn,  ana  the  part  of  thermal  strain; 
1  def j ,  l.e.  .  iJ 
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The  elastic  strain  components  de|j 
components ,  dojj ,  by 


(1) 


are  related  to  the  incremental  stress 
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Dijkl(d,ij"d,kl“MT<kI) 


(2) 


in  which  6  is  the  thermal  expansion  coefficient,  dT,  is  the  change  in 
temperature,  and  6^1  is  Kronecker  delta.  For  isotropic  material,  the  4th 
rank  material  tensor,  la  defined  by 


Dijkl 


E 

1+v 


fiik4jk  +  ( l+v)(  l-2v)  6ij6kl 
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where  E  is  Young’s  modulus  and  v  is  Polsslon's  ratio  of  material. 

On  the  basis  of  the  Von  Mises  yield  criterion  with  temperature-depen- 
dent  yield  strength  of  a  material,  the  yield  surface  can  be  represented  as:: 
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a^j  is  a  censor  representing  Che  total  translation  of  the  center  of  the 
initial  yield  surface,  and  k  is  related  to  the  uniaxial  yield  stress 
<  -  oy(T)//3. 


In  addition  to  the  yield  condition,  a  constitutive  relation  between 
plastic  strain  increments,  stress,  and  stress  increments  is  required  to 
describe  the  Inelastic  behavior  of  a  material.  The  constitutive  relation 
(flow  rule)  used  in  this  Investigation  is  based  on  Drucker’s  postulate  for 
work-hardening  material  [8].  The  flow  rule  is  given  as: 


dV 


df 


do 


ij 


where  dX  is  a  positive  scalar  quantity.  On  the  basis  of  Drucker's  state¬ 
ment,  this  plastic  strain  increment  tensor  must  lie  on  the  outward  normal 
|to  the  yield  surface  at  the  Instantaneous  stress  state. 


Based  upon  Prager's  kinematic-hardening  rule  [4,9]  with  Ziegler's 
modification  [10]  ,  the  Increment  of  translation  of  the  center  of  the  yield 
surface  is  assumed  to  be  directed  along  the  radius  vector  connecting  the 
center  of  the  yielding  surface  to  the  instantaneous  stress  state,  l.e., 

i  I 

do^j  -  (°tj~a1j)du;  dp>0  (7) 

where  dp  can  be  determined  provided  the  stress  point  remains  on  the  tran- 
lated  yield  surface  during  plastic  flow,  i.e.. 


3f 

3o 


dp  * 


kl 


d\L  +  I?11 


(a _ -a  ) 

mn  mn 


3f 

do 


(8) 


mn 


During  plastic  loading,  the  consistency  condition  requires  that 
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+  4?"  - 0 
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;The  plastic  strain  vector  Cdcfj  is  considered  as  the  projection  of  do*, 
(and  thus  of  do^j)  on  the  exterior  normal  to  the  yield  surface,  where  C  is 
a  material  constant.  Hence,  for  small  Incremental  of  stress  and  strain, 
one  can  readily  find  that 
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*•  Therefore,  the  flow  rule  beceuees 
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If  axisymmetric  defoliation  la  considered,  with  reference  to  the 
cylindrical  coordinates  (r,  0,  s)  the  state  Is  defined  by  the  nonvanishing 

stress  components  {do}^  “  <  dor»  dog,  do8,  dtrz  >  and  strain  components 

fdc)1’  “  <  der,  deg,  des,  dyrz  >.  Then,  the  incremental  stress-strain 

relations  are  found  In  the  following  matrix  fora: 

(do)  -  [D]  {de}  -  BdT  {8} 
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METHOD  OF  SOLUTION 

1 

i  1 

tv,  it.  !  The  projectile  and  all  loading  acting  on  the  projectile  are  considered 
as  axisymmetric.  A  cylindrical  coordinate  system  (r,  6,  z)  is  used  in  this 
analysis.  By  the  assumption  of  axial  symmetry,  all  variables  are  indepen- 
rst  dent  of  angle  8,  consequently,  all  derivatives,  with  respect  to  8  vanish. 

*•'  The  displacement  uQ,  and  the  shear  stresses  Tr0  and  t0z  vanish. 

I 

Due  to  the  complexity  of  the  geometry  of  a  projectile  and  nonlinear 
material  behavior,  the  finite  element  method  [11,12]  was  used  to  conduct 
the  stress  analysis.  The  analytical  approach  used  in  this  investigation  is 
the  incremental  loading  technique,  wherein  at  each  step  of  loading  a  new 
stiffness  matrix  is  formulated,  in  terms  of  the  finite  element  model,  and 
solved  for  incremental  deformations,  stresses,  and  strains.  1 
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To  perform  a  finite  element  stress  analysis,  the  cross  section  of  a 
projectile  is  divided  into  a  large  number  of  small  triangular  and  quadri¬ 
lateral  elements  as  shown  in  Figure  1.  Only  a  few  nodal  points  on  the 
boundary  need  be  specified,  the  remaining  nodal  points  are  obtained  from  an 
automatic  mesh  generating  computer  program.  The  grid  was  partitioned  to 
place  a  finer  grid  at  those  areas  of  the  body  that  are  expected  to  undergo 
large  stress. 


PROJECTILE  CONFIGURATION  AND  MATERIAL  PROPERTIES 

♦ 

The  outline  of  a  30  nn  XM  TP  projectile  is  shown  in  Figure  2.  The 
projectile  is  made  of  three  distinct  parts  with  three  different  materials 
as  shown  in  Figure  2  (the  projectile  body,  nose  cap,  and  rotating  band  are 
made  of  1018  steel,  aluminum,  and  gilding  metal,  respectively).  The 
material  properties  are  given  below:  ; 


Poisson's 


Material  Ratio 

Steel  .29 
Aluminum  .33 
Gliding  Metal  .33 


Young's  Yielding 
Modulus  (psl)  Stress  (psl) 

29x10®  88,500 
10x10®  77,000 
17x10®  40,000 


Ultimate  Density 
Stress  (psl)  lb/ln.* 

90,000  .280 
80,000  .101 
48,000  .317 
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►  The  stress-s train  relations  of  those  materials  are  shown  In  Figure  3< 
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Figure  2.  Projectile  Configuration 


Figure  3.  Stress-Strain  Relation 
Projectile  Materials 
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LOADING  CONDITIONS 
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The  loads  considered  in  this  Investigation  simulating  the  environment^ 
in  the  gun  barrel  during  firing*  The  loading  consisted  of  four  types: 

/  1  \  _ 1  1 .  /  A  \  «  ,  _  "  * 


(1)  propellant  gas  pressure,  (2)  setback  force  due  to  acceleration,  (3) 
rotational  velocity  due  to  the  spin  of  projectile,  and  (4)  forced  displace¬ 
ment  of  the  rotating  band*  The  maximum  value  of  these  four  types  of 
loadings  are  given  by,  Pfflax  -  60,000  psi,  -  4.8xl07  in./sec2,  umax  - 

6*335  rad/sec,  and  Am_y  -  -0.008  inches  for  gas  pressure,  acceleration,! 
spin  rate,  and  radial  displacement,  respectively.  i 

The  frictional  shear  forces  between  j 


the  band  and  the  barrel  were 
neglected  in  this  investigation. 
The  loads  were  applied  in  incre¬ 
mental  fashion  with  the  relative 
magnitudes  at  each  load  point 
simulated  the  physical  interde¬ 
pendence  of  the  loads  at  different 
times  during  the  interior  ballistic 
cycle.  The  detailed  loading  history 
is  given  in  Figure  4. 


Figure  4.  Loading  History 


COMPUTATIONAL  RESULTS 


For  each  Incremental  loading,  the  stresses  and  strains  of  each  j 
element  and  the  displacement  of  each  nodal  point  were  calculated.  The 
plastic  flow  has  been  initiated  within  element  186* of  the  gilding  metal  at  ' 
the  end  of  the  first  step  of  loading.  Based  upon  the  Von  Mines  yield  ■ 
criteria,  the  equivalent  stress  for  each  element  was  calculated.  At  the  ! 
end  of  the  first  step  of  loading,  the  maximum  equivalent  stress  (oe)  in  j 
the  projectile  body  is  69,900  psi  which  is  lower  than  yielding  strength  of  ! 
,1018  steel  (element  174),  and  in  the  gilding  metal  portion  is  40,000  psi  • 
(element  186) ,  which  is  equal  to  the  yield  strength  of  the  gilding  i 
material.  j 

The  computation  results  presented  in  this  section  are  based  upon  the  i 
following  three  major  loading  steps  (shown  in  Figure  4): 
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1)  At  the  end  of  the  12th  step  (A-1.0,  p-0.5,  8-0*5,  and  w2-0.036) 
while  the  forced  displacement  reeaches  the  maximum  in  the  rotating  band. 

2)  At  the  end  of  the  33rd  step,  (A-1.0,  p-1.0,  a-1.0,  and  u2-0.l4) 
while  both  propellant  gas  pressure  and  setback  force  reach  the  peak. 

■  3)  At  the  end  of  the  63rd  step  (A-1.0,  p-0,  a-0,  and  u2  -  1.0), 

while  the  spin  rate  reaches  the  maximum  at  the  muzzle  end  of  the  tube. 

j  The  plastic  range  in  the  projectile  under  the  above  three  loading  ! 
conditions  is  shown  in  Figure  5,  6,  and  7,  as  indicated  by  the  shaded  ; 
area. 
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,  The  maximum  plastic  strains  encountered 

g  in  the  gilding  metal  (element  186)  is 

*£>0*  0.027  in. /in.,  and  in  the  projectile  is 

£  o  0.025  in. /in.  (element  163)  at  the  end 

^  y  of  the  12th  loading  step,  while  the  max- 

^  <  imum  radial  displacement  took  place  in 

c  ..  the  rotating  band.  The  maximum  plastic 

^  g0  strains  were  increased  to  0.0294  in. /in. 

(element  186)  in  the  gilding  metal  and 

^  £  7  _ t0  0.0254  in. /in.  (element  163)  in  the 

0  g^g _ projectile  body,  while  both  gas  pressure 

^  °  ^  _ and  setback  force  reached  the  peak. 

men  When  projectile  reached  the  muzzle  end 

£  0  the  maximum  plastic  strains  are  0.0294 

g  q  _ in. /in.  and  0.0258  in. /in.  encountered 

-rC  'J  «  ....  ..in  the  elements  186  and  163,  respective- 
v  _  ly.  During  the  entire  loading  process 

ij  w  °  the  change  of  maximum  plastic  strain  is 

°  small.  The  location  of  plastic  zones 

***  **  ^  under  different  loading  conditions  is 

important  to  the  projectile  design. 


The  equivalent  stress  along  the 
length  of  the  inner  surface  elements 
of  the  projectile  is  shown  in  Figure 
8.  The  maximum  equivalent  stress 
reaches  88,800  psl,  which  is  greater 
than  the  yielding  strength  (Oy  * 
88,500  psl)  and  smaller  than  the 
ultimate  strength  (ou  -  90,000  psi) 
of  the  projectile  material.  If  the 
investigation  was  based  on  elastic 
analysis,  the  maximum  equivalent 
.'stress  could  reach  a  value  which 
would  be  much  higher  than  the  ulti¬ 
mate  strength  of  the  projectile 
material.  This  is  the  reason  that 
our  Inelastic  analysis  technique  be 
developed,  since  a  more  rigorous 
stress  analysis  can  only  be  obtained 
by  using  nonlinear  theories  of  plas- 
, tlclty  as  developed  in  this 
investigation. 


SEE  NEXT  PAGE  FOR  ENLARGED 
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Figure  8.  Equivalent  Stress  in  the 
Inner  Elements  of  the  Projectile. 
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The  equivalent  stress  distribu¬ 
tion  In  the  projectile  of  high  stres¬ 
sed  region  at  peak  pressure  and  set¬ 
back  is  shown  In  Figure  9.  It  Is 
clearly  indicated  that  the  critical 
region  is  located  from  the  rotating 
band  to  the  inside  wall  of  projectile 
base. 

*  The  deformed  grids  under  diffe- 

1  rent  loading  conditions  are  shown  in 
Figures  10,  11,  and  12.  The  deforma¬ 
tions  shown  in  these  Figures  are  the 
deformations  enlarged  10  times  in 
magnitude . 
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CONCLUSIONS 

On  the  basis  of  Drucker's  flow  rule 
of  plasticity,  the  Von  Mises  flow 
criterion,  and  the  strain-hardening 
and  compressibility  properties  of  a 
material,  a  more  rational  rigorous 
nonlinear  elastoplastic  analytic 
method  has  been  develped  for  anal¬ 
yzing  and  designing  a  projectile 
subjected  to  actual  complicated 
firing  conditions.  In  constrast  to 
the  traditional  designing  and 
analyzing  techniques  (theories  of 
elasticity),  the  complex  phenomenal 
plastic  flow  in  the  projectile  has 
been  taken  into  account. 


This  is  due  to  the  necessity  of  designing  for  maximum  stress  to  achieve  an 
optimal  design  for  reducing  component  weight.  Both  nonlinear  material 
response  and  geometric  nonlinearity  have  been  taken  into  consideration. 
Nonlinearity  of  material  properties  has  taken  into  account  by  use  theories 
of  platicity.  Geometric  nonlinearly  has  been  considered  by  use  the  finite 
element  technique.  An  incremental  loading  procedure  has  been  used  to 
consider  the  actual  firing  environment  of  a  gun.  The  complete  loading 
history  which  includes  propellant  gas  pressure,  setback  force,  and  spin 
rate  of  a  projectile  was  defined  as  function  of  time. 

The  critical  region  is  located  in  the  region  between  rotating  band 
and  projectile  base.  In  this  region  the  equivalent  stress  in  general  is 
above  the  yield  strength  of  material,  however,  it  is  below  ultimate 
strength  of  the  material.  The  plastic  strain  has  taken  place  in  the 
relgon,  however,  maximum  equivalent  plastic-strain  at  inner  surface  is 
below  0.03  in. /in.,  which  is  relatively  small. 

Based  upon  our  investigation  the  plastic  zone  in  the  projectile  has 
been  identified  and  accurately  located  for  each  incremental  of  loading 
which  simulating  the  actual  fire  environment  in  a  gun.  The  developed 
technique  will  provide  a  more  rigorous  stress  analysis  tool  for  projectile 
design.  The  potential  design  flaws  and  critical  regions  of  a  projectile 
can  be  Identified  before  the  projectile  is  being  made  and  tested. 


*CHU  and  STEINER 
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